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1. Introduction 

We consider the Wishart-Laguerre (L) and Jacobi (J) ensembles of ran- 
dom matrix theory characterized by the following joint probability density 
(jpd) of real eigenvalues [1, 2, 3]: 

N 

7f\x„ . . . , A^) := Ci,% n 1^. - ^^1' n Aj^^^^^"^e-i^^^ (1) 

j<k j=l 

N 

j<k j=l 

In the above equations, cj^^ ^ and ^ ^ are normalization constants, 
while the index /3 = 1,2, A characterizes the symmetry class of the ensemble 
(orthogonal, unitary and symplectic respectively). 

The Wishart-Laguerre ensemble contains covariance matrices of the form 
W = XXt, where X is a iV x Af (M - TV := i/ > 0) matrix with i.i.d. 
Gaussian entries (real, complex or quaternionic variables). The matrix W 
is symmetric and positive semidefinite, so its N real eigenvalues are non- 
negative ({Aj} > 0). Originally introduced by Wishart [1], matrices from 
this ensemble have been extensively used in multivariate statistical data 
analysis [4, 5] with applications in various fields ranging from meteorological 
data [G] to finance [7, s]. They are also useful when analyzing the capacity 
of channels with multiple antennae and receivers [9], in nuclear physics [10], 
chiral quantum chromodynamics [11] and also in statistical physics such as 
in a class of (1 + l)-dimensional directed polymer problems [12]. Recently, 
they have also appeared in the context of knowledge networks [13] and new 
mathematical results for the case v < have also been lately obtained [14, 
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15]. Large deviation properties of the eigenvalues have been investigated in 
[16, 17, 18], while for an excellent review we refer to [19]. 

The Jacobi ensemble contains combinations of two N x N Wishart- 
Laguerre matrices Wi and 'W2 of the form: 

a = (Wi-W2)(Wi+W2)-^ (3) 

and its eigenvalues are real and lie on the support —1 < {Xj} < 1. Matrices 
distributed according to the Jacobi weight arise as l) truncations of Haar 
orthogonal, unitary or symplectic matrices (for the case of unitary matri- 
ces, an important application arises in the theory of electronic transport in 
mesoscopic systems at low temperatures as detailed in Appendix A); ll) as 
composition of projection matrices [20]. 

In order to proceed, we first define a shifted version of the Jacobi en- 
semble with eigenvalues between and 1: 

N 

Tf\x^, ...,Xn):= Ci^2a,^ n 1^^- - ^'^1'' n ^"(1 - (4) 

j<k j=l 

which appears more frequently in physical applications (see Appendix A) 
and numerical algorithms [21, 22]. By changing variables Xj = 1 — 2yj in 
(2), it is easy to see that: 

N,l3,a,b 2Ar(a+6+l)+§Af(Af-l) ■ ^ ^ 

We also define the average spectral densities {p{x) = (^jf J2^=i ^(^ ~ 
for the three ensembles above as the marginals of their respective jpds: 

(^) i' \. \ _ / J\ J\ .,a)(i)/ 



PNSui>^l)= dX2---dXNTp>{Xi,...,XN) (6) 



PNPabi^l) = / dX2---dXNr^"iXu...,XM) (7) 



PNLbi^l)= / dX2---dXN'?r'iXi,...,XN). (8) 

It follows immediately from the definition that the above densities are all 
normalized to 1, and that the following relation holds between the Jacobi 
and the shifted-Jacobi densities: 

pKU(1-2x) = ^pS^J,,,(x). (9) 

The purpose of this paper is twofold: 
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• to collect and present in a user-friendly way (well-known, but somehow 
scattered throughout the literature) explicit formulae for the above 
densities for finite and all symmetry classes; 

• to use these formulae to compute integer moments of the eigenvalues 
Tn = {Y2iLi -^F) (where the average is taken w.r.t. any of the three 
jpds above). More precisely, we define: 



Ti'\N,p,a,h) 



T(^-^)(iV,/3,a,b) 



One application to the case of quantum transport in chaotic cavities 
is detailed in Appendix A. Other interesting mathematical results for 
the Laguerre case can be found e.g. in [23] and references therein. 
We are also aware that formulae for the Wishart-Laguerre and Jacobi 
moments for /3 = 1,2,4 and finite have been derived by Mezzadri 
and Simm [24]. These are different from the ones provided here and 
obtained via a different method, but equivalent. 

It is easy to see that the average of any linear statistics (i.e. a quantity of 
theformyi = EL/(Ai)), 

{A) = j dXi--- dAjv^MAi, ...An) f{\)^ (13) 

can be computed as a one-dimensional integral over the corresponding spec- 
tral density of the ensemble as: 

N 

(E/(^«))=^ / dxp{x)f{x). (14) 
i=l 

The technical achievement we report in this paper is an explicit computation 
of this integral valid for finite matrix dimension N and all three /3's for the 



/ d\i---d\N'yf\\i....AN) (E^n 

(1 

( dXi---dXN'Pi^\Xl,...,XN) ifZ^i] 



(11) 



[ d\,---dXN7'f\xu...,XN) (y^xA. 



(12) 
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case /(x) = x", n G N. The integral (14) also elucidates a possible strategy 
to evaluate a regular asymptotic expansion of moments for large A^, which 
has been recently highlighted [25] as a problem of current interest in the 
context of electronic transport in chaotic cavities (see Appendix A): all 
one has to do is to seek for a regular (1/A^) expansion of the macroscopic 
spectral density of the form: 



p{x) = p(-)(x) + (x) + -^p(2)(x) + . . . (15) 



in the spirit of high genus correlator expansions [20], and then integrate 
term by term. While this program is far from completion and is thus left 
for future work, we show in Appendix B that at least the leading order 
of the asymptotic expansion of moments is well reproduced for the case of 
quantum transport in cavities with broken time-reversal symmetry (/3 = 2). 

The plan of this paper is as follows: in Section 2 we consider the Wishart- 
Laguerre ensemble. We first summarize the densities for all three 13s, then 
use the integral formula (14) to compute integer moments for /3 = 1,2,4 
in the three subsections. Then in the last subsection we collect results 
from numerical simulations for the density and moments. The same thing 
is done for the Jacobi ensemble in Section 3. A summary and outlook is 
provided in Section 4. In Appendix A we provide a detailed introduction 
to the problem of quantum transport in chaotic cavities which constitutes 
the main motivation for this study, while in Appendix B we show that the 
leading order term in the expansion of the moments via spectral density is 
correctly reproduced. 



2. Wishart-Laguerre ensemble 
2.1. Spectral densities 

The spectral density for the Laguerre ensemble is known for all three 
symmetry classes [27, 28, 29, 30] and, after tedious algebraic manipulations 
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can be cast in the form p'^'^p = i^'J^^'^^ ^{x /2), where: 

m=0 

■^Ar,4,i/V-^/ — 2 2Af,2,2;/V-^i ^ ^A^ 22!^+lr(iV + + 1/2) 

^(r?^ + l/2) (2.) , . , . 



where: 



r(m + i/+ 1)' 

m=0 



(«;+7V-2)/2 

<^i(x) = (2x)^e-2- Yl f-4:i+i-«(2x) (19) 



m=0 

r(m + 1 - k/2) 



2r((. + i)/2,x) ■ 
''^ r((^ + i)/2) ^. 



2'^-ir(m + (z^ - l)/2 + 2 - k/2) ' 



(20) 
(21) 



In the above equations, k = N mod 2, L^^\x) is a generahzed Laguerre 
polynomial defined by the sum: 

n 

LW(z)=5^c,(n,A)z^ (22) 
fc=o 

where: 

r(A + n + l)(-n)fc 
n!fe!r(A + fc + l) 

and r(a, x) = J"^ t°'~^e~^dt is the incomplete Gamma function, while (x)„ = 
r(x + n)/T{x) is the Pochhammer symbol. 

In the following subsections, we present the results for integer moments 
t^\n^ = SjLi -^^^ obtained by integration (eq. (14)) of the densi- 
ties given above. The task is most easily accomplished by first defining the 
following auxiliary function: 

/•oo 

Q(r;m,i;a) := dx x"" e-"" L^^\x)Lf\x) (24) 
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which is easily evaluated using (22) as: 

m I 

Q(r; m, l\ a) 

fc=Ofc'=0 

Note that, by the orthogonality relation of Laguerre polynomials, one has: 



(£ + a)! 
Q{a;m,i;a) = 6mi- 



2.2. Moments (3 = 1 
For /3 = 1, the final result of the integration reads as follows: 



(26) 



)(iV, 1, = 2-"rW(iV, 2, u) + J2(n) + J^in) 



(27) 



where: 



r((iv + i)/2) 

2r((iV + z.)/2) 

_ r((jv + i)/2) ' 

2r((iV + iy)/2) 



{k+N-2)/2 



^ t)„Q(n + i/;iV - 1,2m + 1 - k;i/) 



m=0 



(1-k) 



2{3-i.)/2 



r((z. + i)/2) 

c,„(iV - 1, z/)2"+(^+i)/2+™r(n + z/ + 1 + m 



yi(n)-2(i-^)/2y2(n) 



(28) 



+ 



(29) 



n+{u + l)/2 + m 



N-l 
m=0 

X 2-^1 (ra + (i^ + l)/2 + m,n + z^ + l + m;n + (i/ + 3)/2 + m; -1) 

(30) 

y2(n) = ^ c„(Af - 1, i/)2"+('^+i)/2+mp(^ + m + (z. + l)/2). (31) 

m=0 

In (30), we have used the following hypergeometric function: 

Z7 ^ ?, \ (ai)fe(a2)fc fc 

2Fi(ai,a2;6i;z) -=2^^^^ — —z . 



k=0 



kl 



(32) 



One can check by direct inspection that Tq^' {N, = N and t[^^ {N, 1, v) 
N{N + i^) as it should be. 
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2.3. Moments (3 = 2 
Combining (14) with (17), one easily obtains: 



r(m + 1) 



N-1 



^-^ Tim + v 

m=0 ^ 



1) 



Q{n + v; m, m; v) 



(33) 



One can check by direct inspection that Tq^^ {N, 2,v) = N and t[^'^ {N, 2, u) 
2N{N + I/) as it should be. 

2.4. Moments (3 = 4. 
Similarly for /3 = 4 one gets: 



N-l 2N 2m 



f\N,4,.) = l4^\2N,2,2u) -Y^Y^E 



Xm,k,k' 



m=0k=0k'=0 



where: 



fm{N, u,n,k,k' 



Xm,k,k' = fm{N,u,n,k,k')ck{2N,2v)ck'{2m,2u) 

T{l + k + k' + n + 2u)T{N + l)T{m + 1/2) 



22^+i-nr(iV + V + l/2)r(m + V + 



(34) 

(35) 
(36) 



One can check by direct inspection that Tq^' {N, A,v) = N and r^^^ ( A^, 4, v) = 
AN{N + u) as it should be. 

2.5. Comparison with numerics 

In Fig. 1 we plot the analytical formulae (17), (16) and (18) together 
with numerical diagonalization of matrices W from the Laguerre ensemble 
with /3 = 1,2,4 respectively^, obtained in Matlab as follows^: 

(3 = 1 X=randn(N,M) ; W = X*X' ; 

(3 = 2 X=randn(N,M)+i*randn(N,M) ; W = X*X' ; 

(3 = 4: X = randn(N,M)+i*randn(N,M) ; 

Y = randn(N,M)+i*randn(N,M) ; 
A = [X Y; -conj(Y) conj(X)]; W = A*A' ; 

In Table 1 below, we compare the analytical results for moments of the 



^ Note that for /3 = 4 the matrices corresponding to a certain TV actually have size 
2N X 2N and thus have 2N real and positive eigenvalues. Only of them are distinct 
though, and only those must be used when carrying out numerical simulations for 
moments. 

^ Alternatively, one can use the tridiagonal algorithm by Dumitriu and Edelman [21]. 
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Fig. 1. Average spectral density of the Laguerre ensemble ioi N — 8,1^ = 7 and 
f3 ~ 1,2,4 (top to bottom). Numerical diagonalization is given in red dots, while 
theoretical results are in solid black. 



Laguerre ensemble with numerical simulations via the Matlab algorithm 
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given above. 



Tabic 1. Comparison between tlieory and numerics for moments of the Laguerre 
ensemble {N = 5,M = 7^i^ = 2) (the nmnerieal results are obtained averaging 
over O(IO^) samples). 



/3 


n 


Theory 


Numerics 


1 


1 


35 


34.994 


1 


2 


455 


455.27311 


1 


3 


7665 


7660.026 


2 


1 


70 


69.981 


2 


2 


1680 


1680.238 


2 


3 


50400 


50409.428 


4 


1 


140 


139.982 


4 


2 


6440 


6437.074 


4 


3 


362880 


362857.134 



3. Jacobi ensemble 
3.1. Spectral densities 

The eigenvalue densities for the Jacobi Ensemble (x S 
[29, 30]^: 



-1,1]) read 



N/2-1 



^7V,l,a,6(^) = ^ J2 *^^2m) ^ [(t>2ni{x)'ll^2m+l{x) - (t)2ni+l{x)'ll^2m{x)] (37) 



m=0 



N-1 



(38) 



n=0 



SJ) _ iJ) 



N-1 



2N,2,2a+l,2b+l 



(^\ _ p(2a+l,2;,+l). X v;- ^(a,b) p(2a+l,2fe+l). 
yx) I'2j\f {xj ^2N,m^2m, 



m=0 



(39) 



where: 



^ For /3 = 1, A'^ is even. The expression for A'^ odd is slightly more comphcated ['2!)]. 
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^ 2"+^+^ r(n + a + l)r(n + 5 + l) 

,2a+l,2fe+l 

fi'2m = 52m+l = /i2m (41) 
02„^(X) = Wa,b{x)PS:^''''^'\x) (42) 

V'2m+l(2;) = 'U;a+l,fe+l(a;)-f2m'''^'^^^"^^(^) (^3) 

<^2n^+l(x) = ^.,,(X) [A2™+lP£tl''''^'^(^) " i?2m-lPi^"+i '''^'^ (x)] (44) 
1 

^^2771(2;) = 2 J sign(x - y) (t)2m{y) (45) 

(a,6) _ (4m + 2a + 2b + 3) r((jV + 2)/2) T{{N + 2a + 26 + 4)/2) 
iv,m - 22«+2b+3 r((iV + 2a + 2)/2) r((7V + 26 + 2) /2) ^ 
^ r(m + l/2) r(m + a + 6 + 3/2) 

r(m + a + 3/2) r(m + 6 + 3/2) ^ 



. a-2 

a = 

2 

where sign(z) = z/|2;| and: 



(47) 
(48) 



^ ra(n + 2a + 26 + 2) _ ^ (ra + 2a + 2)(n + 26 + 2) _ _ 
2n + 2a + 26+l ' "~ 2n + 2a + 26 + 5 

(49) 

In the formulas above, Wa^h{^) = (1 ~ x)"(l + x)^ is the Jacobi weight 
function, and P^f"'^^ is the n— th order Jacobi polynomial with parameters 
a and 6, defined as: 

e"(^)^|:("r)(:^')(^)"l^)'- 

1 " 

= ,^,Y.'^f\a,b){x-ir-'{x + iy (50) 



j=0 



where we set: 



Jn). ^ r(n + a+l) T{n + b+l) 

r(i + i)r(n + a - j + 1) r(n - i + i)r(6 + i + 1) • 
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We now turn to the problem of computing integer moments of Jacobi matri- 
ces. It is easy to derive first the following relations between integer moments 
in the shifted- Jacobi and ordinary Jacobi ensembles: 



.{sJ), 



r--(iV,/3,a,b) 



^ \k 

k=0 



2""y: ? (-l)M'^(A^,/3,a,b) 



(52) 



4-^) (iV, /3, a, b) = Q {-2)'4''^ {N, /?, a, b). 



In the following we will therefore focus on the ordinary Jacobi case. 

3.2. Moments (3 = 1 
The n-th moment of the distribution is therefore given by: 



iV/2-l 



2m) 



-1 / r(") 



m=0 



L 



in) 

2m+l,2m 



where: 



An) 



dx (j)j{x) ipki^) x^' 



(53) 



(54) 



(55) 



are computed explicitly in (58) and (63) below. As it is clear from equation 
(54), we only need to compute two different kinds of integrals. Let us then 

start from the first kind {I2m2m+i)- 

Using (42), (43) and (50), we therefore get (in the following the cj""^ 
coefficients will always depend on the pair (2a + 1, 2& -|- 1), so we shall omit 
their explicit dependence on those parameters for the rest of the subsection): 



Tin) 

'2m,2m+l 



^ 2m 

dx (l)2m{x) ^2m+i(a;) 2;" = ^ X] 



Am-i-j i'ini) (2m) , 



' C ■ C ■ 



i,j=0 



-1 



dx (1 — X 



|4r?i+2a— i— jr + 



(56) 



£(")(4m+2a-i-j+l,2fe+i+i+l) 

where we have introduced the following integral: 
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j2i''\-'2fB{w + k + l,z + l) (57) 

I,— n \ / 



X 

fc=0 



introducing Euler's Beta function B(x,y) = V{x)T{y) /T{x + y). Inserting 
this result into equation (56) we obtain: 



_ 92(a+b+l)+l / n4m-i-j J2m) {2m) 

^ -'2m,2m+l " ^ V -LJ C- C- X 



i,i=0 

n 

X 



(t) (-2)''B (2(2m + a + 1) + A; - i - j, 2(6 + 1) + i + j) . 

(58) 



In order to compute the second type of integrals (^2™+! 2m) ™ (^4) ^'^'^ 
(55), we use (44), (45) and (50) to get: 



2m 

j{n) 

-'2m+l,2m 22m+l 

i=0 



;^5^cf") / dyu;,,,(2/)(y-l)2™-(y + l)^x 

r+l 

J dx sign(x - y) (j)2m+i{x) . (59) 



The integral 3^2^+1 (y) '^^'^ ^e rewritten as follows by exploiting (44): 



3^2m+i(2/) = ^2m+i9l';^(2m + 1; y) - i?2m-iSl;i(2m - 1;?/) (60) 



where we have: 
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1 ™ r+l 

= ^ J](-l)™-^cf ^ / dx sign(x - y) (1 - ^ (1 + = 

= 2'^+6+i ^(_i)™-,,M ^ r (-2)'= Bi-, (a + m + - i + 1, 6 + j + 1) 

o— n 1 — n ^ ^ 



jr=0 fc=0 



(61) 



where By{w,z) = 2 By{w, z) -B{w, z) and By(zi;,z) = f^dt t'"-i(l - 
is the incomplete Beta function. Therefore we get: 



2m+l 

^2m+i J] (_i)2m-i+i^(2-+i) /2(m + 1) + a + - j, 6 + j + 1 



X 



2 

j=0 



2m-l 

N|2m-i-l^(2m~l) ^ 

i=o 



- S2n^-l J] (_l)2m-,-l^(2m^l) + a + - j, 6 + j + 1 



(62) 



Eventually, inserting this result into (59) and computing the integral in dy, 
we obtain 



(n) _ 2-+^ (2„^) V l^'' V-2)^ X 

^ -'2m+l,2m - 22m "^j I i 

j=0 fc=0 ^ ^ 

2m+l 

A2m+i Yl i-lf^-^+'cf'"-''^ na,b{2m-i,i;2{m + l) + a + k-j,b + j + : 
j=0 

2m-l 

-B2m-i ^ (-l)2™-i-icf"^-') ^a,b[2m-i,i-2m + a + k-j,h + j + l'^ 

j=0 

(63) 

where we have: 
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/ dx w„ h(x)(x - l)^(x + lY Bi-x {w,z) 



h r)a+b+h+e+l 



B(l + b + i,l + a + h + w)x 



-If 2 



X 3F2(w,l + a + h + w,l - z;l + w,2 + a + b + h + £ + w;l 



B{w,z) B[l + a + h,l + b + - 



(64) 



In (64), we have used the following generalized hypergeometric function: 

(65) 



/ A \ («i)fc(a2)fc(«3)fc k 
iF2[ai,a2,a3;bi,b2;z) ■= — —z . 



k=0 



{bi)k{b2)k k\ 



3.3. Moments (3 = 2 

We wish now to compute the generic n— th moment of the density in 
(38), which reads 



X 



(66) 



By making use of the representation of Jacobi polynomials in (50), the 
n— th moment in (66) can be written as follows (throughout all the present 
subsection the cf^ coefficients will depend on the pair {a,b)): 



T^'\N,2,a,b)=^ih 



fc ^)~^ J dx Wa,b{x) 



k=0 

N-1 k (k) (k) 
k=0 i,j=0 ^ "fc 



(67) 



where: 



&'Hk,i,j,a,b) = (-l)2fc-*-i / dx (1 -x)2*^+"-^-^(l + x)''+*+^x" (68) 
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and can be computed by changing variables and setting y = (1 — x)/2. 
When doing so, one obtains: 



l'^''\k,i,j,a,b) = (_i)2fc-i-i22'=+'^+''+i / dyy2fc+a-i-i(i_y)6+i+i(i_2y)'^, 

Jo 

(69) 

and expanding (1 — 27/)" according to the binomial Theorem, one gets: 

I(-Hk,^,j,a,b) = (-1)2'=— ^22'=+«+^+if](-2)^(';)x 

i=o 

Jo 

The integral appearing in the previous equation is of the following kind 

C dte-^{l-t)y-^ = B{x,y). (71) 
Jo 

Thus we get: 



c—n ^ / 



X 



t=0 

yi^{2k + a + l-i- j + l,b + i + j + 1). (72) 

Plugging this result in equation (67), one can easily see that the n— th 
moment has the compact expression: 



iV-l k n 

ry)(iV, 2, a, 6) = ^ ^ ^ t)fc,^,,/B(2A; + a + £ - i - j + 1, 6 + i + j + 1) 

fc=0 ij=0 1=0 



(73) 

where: 



IP 



= — Thr^ ( ^ ) • (^^) 



After inserting (73) into (52), we have checked that the special case o = 

/i, b = for Tn'^\N, 2, 0, b) indeed agrees with eq. (16) in [31] and with eq. 
(13) in [32] as it should (see Appendix A for details). 
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3.4. Moments (3 = 4: 
The n-th moment of the density in (39) reads: 



r(^) (N, 4, a, b) = ^r^-^) (2iV, 2, 2a + 1, 26 + 1) - i 9 



(75) 



m=0 



where: 



0^,- = y&rni'^^ + 1> 26 + 1) (76) 

= dx Pk''\^) Pt''\^) (77) 

Exploiting the Jacobi polynomial expansion (50) and the binomial Theorem 

we get (again, we omit the explicit dependence of the c- coefficients on the 
pair (a, 6)): 



1=0 j=0 

r+1 

X J dx {1- x)'^+^+^-^-J (1 + x)''+'+^ = 

k £ n y \ 

i=0 j=0 s=0 ^ ^ 
xB{a + k + e + s-i-j + l,b + i + j + 1). (78) 



3.5. Comparison with numerics 

In Fig. 2 we plot the analytical formulae (37), (38) and (39) together 
with numerical diagonalization of matrices from the Jacobi ensemble with 
/3 = 1,2,4 respectively, obtained in Matlab using the algorithm by Edelman 
and Sutton [22]. 

In the following tables, we compare the analytical results for moments 
of the Jacobi and shifted-Jacobi ensembles with numerical simulations via 
the algorithm by Edelman and Sutton [22]. 




Fig. 2. Average spectral density of the Jacobi ensemble for = 6, a = 5, 6 = 9 and 
/3 = 1,2,4 (top to bottom). Numerical diagonalization is given in blue triangles, 
while theoretical results are in solid black. 
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Table 2. Comparison between theory and numerics for moments of the Jacobi 
ensemble {N = 6, a = 5, & = 9) (the numerical results are obtained averaging over 
0(10^) samples). 



(3 


n 


Theory 


Numerics 


1 


1 


1.1428 


1.1454 


1 


2 


1.1536 


1.1573 


1 


3 


0.5108 


0.5125 


2 


1 


0.9230 


0.9259 


2 


2 


1.4872 


1.4892 


2 


3 


0.5470 


0.5485 


4 


1 


0.6666 


0.6674 


4 


2 


1.9300 


1.9294 


4 


3 


0.5602 


0.5602 



Table 3. Comparison between theory and numerics for moments of the shifted- 
Jacobi ensemble {N = 6, a = 5, b = 9) (the numerical results are obtained averaging 
over O(IO^) samples). 

/3 n Theory Numerics 

1. i 2.4286 2.4300 
1 2 1.2170 1.2169 
1 3 0.6902 0.6902 



2 


1 


2.5385 


2.5390 


2 


2 


1.4103 


1.4102 


2 


3 


0.8932 


0.8919 


4 


1 


2.6667 


2.6669 


4 


2 


1.6492 


1.6493 


4 


3 


1.1605 


1.1612 



4. Conclusions 

In conclusions, we have first collected well-known formulae for the spec- 
tral density of Laguerre and Jacobi ensembles of random matrices with or- 
thogonal, unitary and symplectic symmetry. We feel that this paper might 
serve as a quick reference point for formulae that might be hard to dig out 
in the literature. Using these results, we have computed the average of 
integer moments, reducing the complexity from a A^-fold integration to a 
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single integral over the spectral density. In all cases, expressions different 
from ours and derived through different methods already exist [24, 31, 32] 
(see also Appendix A). It would be interesting to prove mathematically 
the equivalence of various formulae which are now available for the same 
objects. The present paper offers a different and possibly simpler way of 
deriving transport moments for all symmetry classes and finite number of 
open electronic channels. The obtained results have been checked numer- 
ically with high accuracy, and the corresponding Matlab codes are freely 
available on demand. 
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A cavity of submicron dimensions etched in a semiconductor can be con- 
nected to the external world by two leads supporting Ni and N2 electronic 
channels. A voltage difference V applied between the two leads lets an 
electronic current flow through the cavity, whose intensity presents time- 
dependent fluctuations which persist down to zero temperature [34]. These 
are associated with the granularity of the electron charge e. Typical phe- 
nomena observed in experiments include weak localization [35], universality 
in conductance fluctuations [30] and constant Fano factor [37]. According 
to the Landauer-Biittiker scattering approach [34, 38, 39], the statistics of 
quantum transport observables depend in a rather simple way on the scat- 
tering matrix S of the cavity. This is a unitary A^o ^ ^0 matrix (where 
A^o = Ni + N2) which relates the wave function coefficients of the incoming 
and outgoing electrons in a natural basis: 



The transmission (t,t') and reflection (r, r') blocks are submatrices encod- 
ing the transmission and reflection coefficients among different channels^. 
The eigenvalues of the hermitian transport matrix T = tt^ are of primary 
importance: for example, the dimensionless conductance and the shot noise 

* (t,t') are respectively of size A^i x A''2 and x Ni, while (r, r') are of size 7V2 x N2 
and iVi X iVi . 



Appendix A 

Electronic transport in open cavities 




(A.l) 
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are given respectively by G = Tr(T) [38] and P = Tr[T(l - T)] [-ID], while 
other transport properties are encoded in higher moments r„ = Tr[T"']. 

When the classical electronic motion inside the cavity can be regarded 
as chaotic, Random Matrix Theory has been very successful in describing 
the statistics of universal fluctuations in such systems: the scattering matrix 
S is assumed to be drawn at random from a suitable ensemble of matrices, 
with the overall constraint of unitarity [41, 42, 43]. A maximum entropy 
approach (under the assumption of ballistic point contacts [34] ) forces the 
probability distribution of § to be uniform within the unitary group, i.e. S 
belongs to one of Dyson's Circular Ensembles [27, 44]. 

The uniformity of § within the unitary group induces the following 
remarkably simple jpd of transmission eigenvalues {Tj} of the matrix T 
[34, 43, 45]: 

'?f\T,, . . . ,T;v) « n l^i - Tkl^W^^^^'^'^"' (^-2) 

j<k 1=1 

with N = mm{Ni, N2), fi = \Ni — and /3 = 1,2,4 in case of preserved 
time-reversal symmetry, broken time-reversal symmetry and spin-flip sym- 
metry respectively. The jpd (A. 2) is precisely of the shifted-Jacobi form (4) 
with a = f (/i -Fl) - 1 and b = 0. 

The transmission eigenvalues are thus correlated real random vari- 
ables between and 1, whose jpd (A. 2) in principle allows for a complete 
characterization of statistical properties of experimental observables. For 
most recent analytical results, we refer to [30, 31, 32, 33, 46, 47, 48, 49, 50, 
51, 52, 53]. 

In particular, the study of higher moments of the transmission matrix 
r* = Tr[T"'] has recently seen many analytical progresses^ [24, 25, 30, 31, 
32, 33, 48, 49, 53]. In particular, we now have two different (but equivalent) 
formulae for higher moments for f3 = 2 and arbitrary A'^i , '■ 

N-l p 
p=0 k,l=0 ^ 



9p{k)gp{t) 
n + k + £- 

n = E t}l(--') see [32] (A.4) 



n-l 
p=0 



where: 



In our notation, r* = r^^'''(A*', 2, /i, 0). 
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Novaes [33] has also computed the leading O(A^) term in the large Ni,N2 ^ 
1 expansion of the moments Tn'^\oo,2,iJ,,0) (Eq. (A. 4)) for the quantum 
transport problem at /? = 2. His formula in our notation reads: 



r(^-^)(oo,2,/i,0)~(iVi + iV2)j; 

p=i 



1 



(A.6) 



where Cp = ^(jf) and ^ = NiN2/{Ni + N2)'^. In the next Appendix, we 
will show how this result (for the case /i ~ 0(1) = 1/4) can be derived 
by integration of the asymptotic spectral density of shifted Jacobi ensemble 
with a = b ~ 0(1) when — )• oo. 

In Fig. 3 we plot the transport moments Tn{Ni,f3, {(3/2){fj. + 1) — 1,0) 
for n = 3 and n = 4 as a function of Ni, while the number of channels 
N2 = Ai + 2 is held fixed. All cases yield expressions formally different 
from others [24] (see equation (52) and the equations in subsections 3.2 and 
3.4), but equivalent. 



o 



4 
3.5 

3 
2.5 

2 
1.5 

1 

' 0.5 





,1,1/2,0) 




,1,1/2,0) 


-0 


,2,2,0) 




,2,2,0) 




,4,5,0) 




,4,5,0) 




1 



10 



12 



Fig. 3. Transport moments r„(A^i, /?, (/3/2)(/i + 1) — 1,0) for n = 3 (circles) and 
n = 4 (triangles) as a function of Ni . The number of channels iV2 = iVi + 2 is held 
fixed. The results are for /3 = 1 (short-dashed red), /3 = 2 (long-dashed blue) and 
13 = A (solid black). 



Appendix B 

Large N asymptotics for moments 

The general integral formula (14) offers a neat and explicit way to obtain 
the large N asymptotic behavior for moments and possibly the full 1/A 
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expansion as follows. Suppose that the spectral density of the ensemble 
under discussion admits the following series expansion: 

p{x) = + lp(i)(x) + ^p^'\x) + ... (B.l) 

then the corresponding series for moments can be formally obtained as: 

Tn = N j dxp^^\x)x'' + j dxp^^\x)x'^ + ^ j dxp^^\x)x'' + ... (B.2) 

A first check of the validity of the putative formula (B.2) goes as follows. 
Novaes [.33] has computed the leading 0{N) term in the large Ni,N2 S> 1 

expansion of the moments Tn '^\oo, 2, p,, 0) for the quantum transport prob- 
lem in chaotic cavities with unitary simmetry supporting Ni,N2 electronic 
channels in the leads {p = \Ni — and = min(Ai, N2)), see eq. (A. 6). 
Now, the A — )• 00 limit for the spectral density of the ordinary Jacobi 
ensemble with a, 6 ~ 0(1) reads [29]: 

vrv 1 — x'^ 

i.e. it is independent of a, 6, /3. Using (9) we get: 

and using (B.2), the 0(A) term in the expansion of the n-th moment of the 
/3 = 2 shifted-Jacobi ensemble should read: 



hm lr(^^)(A,2,a,b)= 



2x'- 

dx- 



TTyJl - (1 - 2xf 

r(n + l/2) 



\AFr(l+n) 



(B.5) 



Eq. (B.5) should then be compared with (A. 6) for a special value of ^, 
namely the one that guarantees that parameters a, b are finite (i.e. of 0(1)) 
when A — )• 00 (otherwise eq. (B.3) would not hold). Since a = ^ = N1 — N2, 
we need to impose that this difference is finite when Ai,A2 — )• 00. This 
immediately leads to ^ = 1/4 and Ai + A2 = 2A. Inserting these values 
into (A. 6), we indeed find the identity: 

^vp-iy 0rr(l + n) 
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Therefore, we have proven that a simple one-dimensional integral over the 
asymptotic spectral density of the shifted-Jacobi ensemble correctly repro- 
duces Novaes' formula for the asymptotics of transport moment (obtained 
by a totally different method) when the difference between the number of 
electronic channels in the two leads is of 0(1). 
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